Spin diffusion of lattice fermions in one dimension 
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We study long-time spin difTusion of harmonically trapped lattice fermions in one dimension. 
Combining thermodynamic Bethe ansatz approach and local density approximation, we calculate 
spin current and spin diffusion coefflcient driven by the population imbalance. We find spin current 
is driven by susceptibility efi^ects rather than typical diffusion where magnetization would transport 
from regions of high magnetization to low. As expected, spin transport is zero through insulating 
regions and are only present in the metallic regions. In the weak coupling limit, the local spin 
diffusion coefficient shows maxima at all the insulating regions. Further, we estimate damping rate 
of difi'usion modes in the weak coupling limit within the lower metallic portion of the cloud. 



I. INTRODUCTION 

The transport of fermions is central to many areas in 
physics. Fermion transport can be seen in various sys- 
tems, including electronic matter, stars during collapses 
and explosions, and quark-gluon plasma. In particular, 
the transport of spin is highly relevant to potential tech- 
nological applications known as spintronics. The subject 
of spintronics is to design devices to control spin evolu- 
tion and dynamics for information storage, transfer, and 
manipulation [Ij] . Thanks to the current state of cold gas 
experiments, cold gas systems provide a unique opportu- 
nity to study these transport properties in a controlled 
fashion. 

The ability of engineering optical lattices in arbitrary 
geometries in different dimensionalities opens up excit- 
ing avenues to visualize and gain deeper understanding 
of fundamental many body physics. Recent progress in 
experimental techniques of ultra-cold atoms, such as the 
technique of detecting a single-site [3, [s^ , allows one to 
probe not only the static properties, but also the dynam- 
ics of the atoms in optical lattices. 

On shorter timescales, a systems of ultracold atoms can 
be used as an effective simulator for thermal equilibrium 
physics in closed many-body quantum systems. In con- 
trast, on longer time scales, the presence of underlying 
inhomogeneous trapping potential provides an ideal plat- 
form for studying numerous non-equilibrium phenomena. 
In cold gas experiments in the presence of population im- 
balance, inhomogeneous trapping potential creates a den- 
sity imbalance throughout the lattice. This density im- 
balance mimics a magnetization and drives the spin cur- 
rent in the cold atom setups. In solid state devices spin 
current is produced by injecting spins into the device and 
can thus be probed by magnetic microscopy or neutron 
scattering. After the recent proposal on the ultra-cold 
atom analog of electronic semiconductor devices known 
as " atomtronics" Q , the study of transport phenom- 
ena gained a tremendous momentum [T-I^, 18|. In the 
absence of optical lattice, spin transport in strongly inter- 
acting Fermi gases has been experimentally investigated 



in two laboratories at MIT and Rice University [^, [13] ■ 
In general, the non-equilibrium dynamics of cold atoms 
in optical lattices are studied after an adiabatic or sud- 
den change of the atom or lattice parameters 12, 19 26]. 
While the dynamics for the adiabatic process are ex- 
pected to follow the change of the system Hamiltonian, 
the sudden quench leads to the non-equilibrium physics 
that allows for the understanding of relaxation dynam- 
ics in the presence of many-body interactions. Although 
non-equilibrium dynamics due to the adiabatic or sud- 
den changes have been widely studied, longer time scale 
dynamics due to the inhomogeneous trapping potential 
have been largely unexplored. 

The spin of the atoms referred to in this paper are one 
of two hyperfine states available to the atoms. There- 
fore, we treat the spin as a scalar quantity which points 
either parallel (up) or antiparallel (down) to some arbi- 
trary direction of quantization. The magnetization M 
defined below is proportional to the difference in up and 
down atom densities. As the overall magnetization in 
the system is produced by the fixed population imbal- 
ance and the interactions are spin- independent, the total 
magnetization is a conserved quantity. The local mag- 
netization cannot disappear locally, but can only relax 
slowly over the entire system. The purpose of this paper 
is to investigate how atoms relax by the magnetization 
being physically transported within the trap. Indeed, 
this transport process occurs only in the presence of in- 
homogeneous potential via a diffusion motion in a larger 
time scale provided by the underlying harmonic poten- 
tial. 

In the present paper, we study the spin diffusion 
of one dimensional lattice fermions in the presence of 
an external harmonic potential. We consider longer 
timescale (defined later) dynamics of population imbal- 
anced fermions and study the spin diffusion current and 
local spin diffusion coefficient, and how these values ef- 
fect the damping rate of the spin diffusion modes. We 
find evidence of spin current in the metallic regions of the 
lattice, while there being no spin current in insulating re- 
gions. Further, in the weak coupling limit, we calculate 
local spin diffusion coefficient for all regions of the lattice. 
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The damping rate of spin diffusion modes are estimated 
using the continuity equation. 

The structure of the paper is as follows. In section II, 
we introduce our model and formalism. In section III, 
we calculate the ratio of spin current to local diffusion 
coefficient for characteristic values of large and small in- 
teraction, and show the interplay between magnetization 
and susceptibility effects to spin current. Further, the 
local spin diffusion coefficient is found over the lattice 
for the weak coupling limit. Finally, the damping rate 
of diffusion modes is estimated within a metallic portion 
of the cloud. We display our results with figures in the 
appropriate dimensionless quantities. In section IV, we 
summarize the results and connect them to experiments. 



II. FORMALISM 

We consider a cloud of atoms consisting of two hyper- 
fine states denoted by spin t and i. We assume that 
the population imbalance P = {N-f — N D/N is finite, 
where N = N-f + N^^ is the total number of atoms in two 
hyperfine states of the same atom. These atoms are sub- 
jected to a combined optical lattice and harmonic trap- 
ping potential. The Hamiltonian of the system can be 
represented by the one-dimensional Hubbard model. 



<ij>,(T i 



(1) 



where cl^icia) creates (destroys) a Fermi atom with pseu- 
dospin a =1,4 at lattice site i. The density opera- 
tor is Ilia- = cl^Cia and < ij > indicates the nearest 
neighbor pair of sites. The average chemical potential 
fi = (/if + lJ-i)/'2 and the chemical potential difference 
h — {^f — iJ,i)/2, where is the chemical potential of 
hyperfine state a. The on-site interaction U can be re- 
pulsive or attractive. Its magnitude and sign can be con- 
trolled by the s-wave scattering length of the two Fermi 
species. In the present work, we consider only the pos- 
itive U. In cold-gas experiments, the ratio U /t is con- 
trolled by the intensity of the standing laser waves (the 
tunneling amplitude t is exponentially sensitive to the 
laser intensity while the on-site interaction U is weakly 
sensitive). Experimentally, the ID geometry can be real- 
ized by a strong confinement in the transverse direction 
with an additional periodic potential applied along the 
other direction. We consider a tight one dimensional ge- 
ometry such that the level spacing in transverse direction 
is much larger than the energy per particle of the axial 
direction. 

In the presence of the underlying harmonic oscillator 
potential, the average chemical potential ^{z) monoton- 
ically decreases from the center to the edge of the trap. 



The spin diffusion current driven by this spatial variation 
of the chemical potential is given by the modified Fick's 
law fii; 
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dz 
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where the magnetization density M(z) — n-f{z) — n^(z), 
the spin susceptibility x{z) = dAI/dh, and the spin diffu- 
sion coefficient D{z) are values that vary in space, even at 
equilibrium. Here n„{z) is the density of hyperfine state 
a at position z. First, we solve the homogeneous Hub- 
bard model for the population imbalance system using 
the thermodynamic bethe ansatz (TEA) technique and 
numerically calculate the magnetization and the spin sus- 
ceptibility as discussed in Refs. [3,[l3l- By combining 
the TEA solutions with the local density approximation 
(LDA), we then extract the local quantities m{z) and 
x{z)- In LDA, the external trapping potential, which is 
independent of hyperfine state, is Vi — at site 

i. This is related to the local chemical potential through 
the relation /ij = /ig — 1/^, where /io is the central chem- 
ical potential and z = id, with lattice constant d is the 
spatial coordinate. We notice that while /i^ is spatially 
dependent due to LDA, h is not and remains constant. 



III. THE RESULTS 

The spin current given in Eq. ([2|) is driven by the den- 
sity imbalance and is a longer timescale result. There 
are several time scales associated with the experimental 
setup. The equilibrium time scale is the maximum oih/U 
or h/t. There are two other larger time scales, one asso- 
ciated with trapping frequency tt/uj and the other asso- 
ciated with the lattice X/{2vp), where is the fermi ve- 
locity. The time scale we consider here is T > tt/lu ^ h/t 
. Systems of trapped fermi atoms without the presence 
of an optical lattice have been studied experimentally 
in Refs. |9|, ll^l, where T > tt/uj, thus, the addition 
of a lattice with the timescale parameter regime being 
T > tt/lu ^ X/{2vp) ^ h/t is feasible in current setups. 



The spin current 

Using LDA, the derivative with respect to spatial co- 
ordinate can be converted into a derivative over the av- 
erage chemical potential /i. Noticing that h is indepen- 
dent of the spatial coordinate and using the relations 
d/dz = {dfi/dz)d/dii and {dji/dz) = —mco'^z, spatial 
dependence of the ratio of spin current and diffusion co- 
efficient is written as 



J,n/D 



(3) 
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Figure 1 shows the results for relatively low temperature 
and large on-site interactions where Mott insulating state 
is present in the middle of the trap. We find that all the 
local quantities, including spin current, show a consid- 
erable spatial variation even in equilibrium (see FIG. [T] 
caption for the details). 

As seen in FIG.[T1 the magnetization is smallest at the 
edge and the center of the trap. The spin current flows 
from low magnetization to high magnetization regions. 
This current pattern originates from two contributions 
as shown in FIG. [2] and the atoms diffuse towards the 
Mott insulating region. On the other hand, the local dif- 
fusion arises from the collisions with opposite spin atoms. 
As the time between scattering event r is proportional to 
the local spin diffusion coefficient (D), D becomes very 
large when the density n and holes (1 — n) are small. 
This is the reason why the ratio of spin current to diffu- 
sion coefficient (jm/D) is almost zero at the center and 
at the edges. The magnetization is almost constant in 
the Mott insulator region, so that (jm/D) is zero in this 
region. A similar current pattern (from edge towards the 
center) has been observed for an attractively interacting 
one dimensional non-lattice fermions 

The results for a weak coupling limit is shown in 
FIG.[3l We find few qualitative differences other than the 
disappearance of Mott insulating region. As one expects, 
Mott region disappears at high temperatures and small 
interactions but the current pattern remains the same. 
We find that the spatial variations of these quantities are 
greater for larger population imbalance (or larger values 
of h). 



The local spin diffusion coefficient 

The fluctuation-dissipation theorem 



x"ik,co)^^{l-e^'^nS{k,uj) 



(4) 



where the inverse temperature /3 = {kBT)~^, and the 
imaginary part of the dynamic magnetization response 
function 



X"{k,uj) = -y dtj cize^-*-''=^([M(z,i),M(0,0)]),(5) 
is related to the magnetization correlation function 



(6) 



Here S{z, t) = {M (z, t)M{0, 0)) and [A, B] represents the 
commutator between the operators A and B. {M) = 



M represents the expectation value with respect to the 
Hamiltonian H . 

First, we use general hydrodynamics description to 
derive the magnetization 27|. Inserting the spin cur- 
rent jVn into the continuity equation dtM + dzjm — 
0, the diffusion equation has the form dtM{z,t) = 
DxV^[M{z,t)/x]- Muhiplying this by e'"* and integrat- 
ing the left hand side, and taking Laplace transformation 
and Fourier transformation, the solution of the diffusion 
equation can be written as 



M(fc,cj) 



iDk' 



:M{k,t^ 0). 



(7) 



Notice that the diffusion process is reflected by the diffu- 
sion pole on the negative imaginary axis a.t uj — —iDk^ . 
As the diffusion coefficient is local, the diffusion pole may 
not be as simple as it looks. 

In order to establish the connection between the hydro- 
dynamic diffusion equation and the dynamic magnetiza- 
tion response function, we assume that spatially varying 
chemical potential difference h{z,t) = h{z)e'^* is adia- 
batically turned on to mechanically induce a non-zero 
magnetization. At time i = 0, the chemical potential 
difference is switched off and allow the induced magneti- 
zation to be relaxed as the system returns to equilibrium. 
The induced magnetization at t = is 



M(z, t = 0) = 2i J dr J dz'x"{z - z', T)e-'"'/i(z'). (8) 

The Fourier transform of this equation has the form 
M{k,t = 0) = x{k)h{k), where 



Xik)= I —X"ik,u;). 

TTUJ 



(9) 



For t > 0, the Laplace transformation of the induced 
magnetization 

M{z,t)=2i I dr I dz'x"{z~ z\t-T)e''''h{z'),{lQ) 

J —OO J 

has the form 

M(fc,.)=/^^4^^Mfc). (11) 

Using h{k) = M{k,t^ 0)/x{k) and 



Xik,v) = 



duj x"{k, 
n (lu — 1]) ' 



we find 



M(fc,7?) = - 



X{k,v) 
Xik) 



M{k,t = 0). 



(12) 



(13) 
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FIG. 1: For strong coupling, the spatial variations of jm{z)/ D{z) (where jm{z) is the spin current and D{z) is the 
spin diffusion coefficient), magnetization M{z) — n-f(z) — n^{z), polarization P{z) = [n^{z) — nj,(z)]/[n^(z) + nj_(z)], 
and atom density n(z) = n-f{z) + n^{z) are shown in panel (a)-(d) respectively. We define the scaled length 
z — z\JmuP' jl. Quantities plotted are dimensionless. We fixed the on-site interaction {U — 5t), the inverse 
temperature (/? = 5/i), and the chemical potential difference {h = OAt). 



Comparing this with Eq. ([7]), the dynamical magnetiza- 
tion response function can be extracted 



iDk^ 



■q- 



iDk^ 



X{k). 



(14) 



Setting rj = Lu + ie and taking the imaginary part at the 
limit of e — ?• 0, we get 



x"(fc,^) = 



-X{k). 



(15) 



Since x"(fc,CLi) is related to the magnetization corre- 
lation function S{k,uj) through the Eq. Q, we no- 
tice that S{k,uj) gives a quasielastic peak with width 
T{k) = 2Dfc2. 

On the other hand, the spin conductivity D{io) can be 
written in terms of the spin current Is{k, t) 



b{uj) = ^ lim -^x"(fc,w) 

X k^o k^ 



(16) 



= — /TOn(w). 

Here the spin current response function is 

U{oj) = --\im dte''^mi,{k,t),U-k,0)]}, (17) 
N k^ajQ 

where the spin current operator Iik,t) = 
X]fc',o- f^sin(fc'd)4,^^Cfc'+fc,cr. Combining Eqs. ([15]) 
and ([TT]). the spin diffusion coefficient can be related to 
the spin conductivity as D = D(0)/2 [2^. Following the 
reference (29j . the long-time spin diffusion coefficient of 
the one dimensional Hubbard model can be casted into 



D 



hxkeT 2 



27rCo 



(18) 



where d is the lattice constant and the coefficients Co = 
(/2(0,0)) and C2 = (/(0,0)[[/(0,0),7J],7J]) are the first 
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FIG. 2: [color online] The dimensionless spin current for 

the same parameters as in FIG. [1] The total spin 
current (blue dashed line) j„i = ji + .12 are sum of two 
contributions. The magnetic contribution (gray line) is 
ji/2D = —dM/dz and the susceptibility contribution 
(black line) is ji/2D = {M/x)dx/dz 



two short-time expansion coefficients of the current corre- 
lation function. Note that the long-time spin diffusion co- 
efficient is written in terms of two lowest order short-time 
expansion coefficients of the current correlation function 

In the weak coupling limit, by generalizing the results 
in Ref. [H, we find 



^ ^ am{kd){nk'\- - n^j.) 



(19) 



y^sin^(fc(i)nfccr(l 

kc7 



nka) 



and 



C2 = U^Cq 



k 



(20) 



Here the Fermi function Uka — {e^'^'" + 1)~^, and eua ~ 
2tcos{kd) — fia- The spatial dependance of the diffu- 
sion coefficient enters through the local chemical poten- 
tial fj.icy = fJ-Oa - 7*^- 

As shown in Fig. 4, we find the diffusion coefficient to 
be small through the middle of the trap, though there is 
a local maximum at the center of the trap as predicted 
since the trap center is occupied by both up and down 
spin particles. A maximum is also seen when the local 
density is equal to unity, a remnant of the Mott transi- 
tion, as well as at the edge of the trap, where the density 



is almost zero, showing that D is large when density and 
holes are small. Further, the spin diffusion coefficient 
tends to reach a local maximum when the magnetization 
is a local minimum, and visa versa. By combining FIG. 
3(a) and FIG. 4, one can see that the spin current is large 
at the edges of the trap in the lower metallic region, finite 
but small in the upper metallic region, and zero in the 
insulating regions. 



The Damping rate 

Taking the inverse Laplace transformation of Eq. Q , 
one can write 



M(k,t) = e-^'' ^M{k,t = 0). 



(21) 



This shows that the damping rate of hydrodynamic dif- 
fusion mode r(fc) = l/(Dk^) is consistent with the 
quasielastic peak of S{k,uj). As the diffusion coefficient 
D = D{k) for trapped fermions, Eq. (PT|) does not give 
us much information on the damping rate of the diffusion 
mode. However, the damping rate of the diffusion modes 
can be estimated directly from the continuity equation by 
assuming M{z,t) has the form M{z,t) = exp[— Tt]M(z). 
Inserting this into the continuity equation and integrat- 
ing both sides, we find 



jm(zf) - j„i(Zi) 

r^M{z)dz • 



(22) 



When estimating r using this equation, one must con- 
sider only a part of the metallic cloud with edges Zi and 
since the spin does not transport through an insu- 



lating region. We estimate r from a small region in the 
lower metallic band, from Zi = 2.21 to zj = 2.42, and 
find T « 2\h/{kBd'^). 



IV. DISCUSSION AND SUMMARY 

The local spin current patterns predicted above for 
neutral atomic fermions in an optical lattice can exper- 
imentally be probed by measuring lattice scale modula- 
tions of the atom density. 



Noise correlations [30|, |3l|, 
Bragg scattering [32l |. and in situ imaging in the lattice 
scaling (ssj are commonly used density mapping tools in 
cold atom experiments. While Bragg scattering is simi- 
lar to X-ray scattering probes of the crystal structures of 
solids, in situ imaging is analogous to scanning tunneling 
microscopy. Imaging the subsequent lattice scale de nsity 
modulations as done for one dimensional Bose gases 3J] , 
local spin current patterns can be probed. 

In this paper we have found evidence of longer 
timescale spin current in both the strong and weak cou- 
pling cases. By calculating the ratio jm{z)/ D{z), we find 
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FIG. 3: For weak coupling, the spatial variations of jm{z)/ D{z) (where jm{z) is the spin current and D{z) is the 
spin diffusion coefficient), magnetization M{z) — n-|-(2;) — n:^{z), polarization P{z) = [n-f (z) — ni{z)\/[n^{z) + n]^{z)], 
and atom density n{z) — n^{z) + are shown in panel (a)-(d) respectively. We define the scaled length 

z = Z\J rrnip' 12. Quantities plotted arc dimensionless. We fixed the on-site interaction [IJ = 0.8t), the inverse 
temperature (/3 = 5/t), and the chemical potential difference {h = 0.2^). 



that the spin current is driven primarily by susceptibility 
rather than directly from magnetization density. In the 
weak coupling limit, we further found how the local spin 
diffusion coefficient varies within the trap, having maxi- 
mums at the insulating rcigions. We find spin current to 
be primarily located edges of the trap. Finally, we esti- 
mated the damping rate of diffusion modes in a metallic 
portion of the cloud. 
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